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The skeleton loop integrals which contribute into the gluon self-energy have been iterated (skeleton 
loops expansion) within the Schwinger-Dyson equation for the full gluon propagator. No any trunca- 
tions/approximations as well as no special gauge choice have been made. It is explicitly shown that 
such obtained general iteration solution for the full gluon propagator can be exactly and uniquely 
decomposed as a sum of the two principally different terms. The first term is the Laurent expansion 
in integer powers of severe (i.e., more singular than 1/g^) infrared singularities accompanied by the 
corresponding powers of the mass gap and multiplied by the corresponding residues. The second 
(perturbative) term is always as much singular as and otherwise remaining undetermined. We 
have explicitly demonstrated that the mass gap is hidden in the above-mentioned skeleton loop 
integrals due to the nonlinear interaction of massless gluon modes. It shows explicitly up when the 
gluon momentum goes to zero. The appropriate regularization scheme has been applied in order to 
make a gauge-invariant existence of the mass gap perfectly clear. Moreover, it survives an infinite 
series summation of the relevant skeleton loop contributions into the gluon self-energy. The physical 
meaning of the mass gap is to be responsible for the large scale structure of the true QCD vacuum. 

PACS numbers: 11.15.Tk, 12.38.Lg 



I. INTRODUCTION 

The Lagrangian of QCD P, |2| does not contain explicitly any of the mass scale parameters which could have a 
physical meaning even after the corresponding renormalization program is performed. However, if QCD itself is a 
confining theory then a characteristic scale has to exist. It should be directly responsible for the large scale structure 
of the true QCD vacuum in the same^ay as Age d is responsible for the nontrivial perturbative dynamics there (scale 
violation, asymptotic freedom (AF) T, 2]). On one hand, the color confinement problem is not yet solved 0. On the 
other hand, today there is no doubt that color confinement is closely related to the above-mentioned structure of the 
true QCD ground state 0i Q and vice- versa. The perturbation theory (PT) technics fail to investigate them. 

The main goal of this paper is to show how the above-mentioned characteristic scale (the mass gap in what follows, 
for simplicity) responsible for the nonperturbative (NP) dynamics in the infrared (IR) region may explicitly appear in 
QCD. This especially becomes imperative after Jaffe and Witten have formulated their theorem "Yang-Mills Existence 
And Mass Gap" @. Moreover, we will show that the mass gap may not only appear, but it may also survive after 
an infinite series summation of the relevant skeleton loop contributions into the gluon self-energy. Thus the paper is 
devoted to a possible solution of one of the important problems in theoretical particle/nuclear physics, namely to the 
dynamical generation of a mass gap in quantum field gauge theories. 

The propagation of gluons is one of the main dynamical effects in the true QCD vacuum. The gluon Green's 
function is (Euclidean signature here and everywhere below) 



D^M = * {TMdii^^O + ^L^M) ^, (1-1) 

where ^ is the gauge-fixing parameter and T^^ici) = 6^^^, — q^q^/q^ ~ S^^ — L^i,{q). Evidently, T^y[q) is the transverse 
("physical") component of the full gluon propagator, while Lf^^{q) is its longitudinal (unphysical) one. The free gluon 
propagator is obtained by setting simply the full gluon form factor to d(Q^,^) = 1 in Eq. (1.1), i.e.. 



D%{q) = * {T^Aq) + ^L^Aq)} ^- (i-2) 

The main tool of our investigation is the Schwinger-Dyson (SD) equation of motion 0,0,13 (and references therein) 
for the full gluon propagator (1.1), since its solutions refiect the quantum-dynamical structure of the true QCD ground 
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state. Some results of the present investigation have been already briefly announced in our preliminary publications 
0, Il0| . Here we present the full (without any estimates, approximations/truncations, specific gauge choice, etc.) 
investigation of the problem of the dynamical origin of a mass gap in quantum Yang-Mills (YM) theory. 

II. GLUON SD EQUATION 

The general structure of the gluon SD equation 0, Q, can be written down symbolically as follows (for our 
purposes it is more convenient to consider the SD equation for the full gluon propagator and not for its inverse): 

D{q) = D\q) + D\q)Tg[D]{q)D{q). (2.1) 

Here and in some places below, we omit the dependence on the Dirac indices, for simplicity, as well as the quark and 
ghost skeleton loop contributions into the gluon self-energy (as it is required by the YM character of our consideration) . 
The nonlinear (NL) pure gluon contribution Tg[D\{q) into the gluon self-energy is the sum of the four topologically 
independent skeleton loop integrals, namely 

UDM = \Tt + \Ti{q) + \T2{q) + \n{q), (2.2) 
where the so-called constant tadpole term is 

Tt^a'j "-§^.nD{q,l (2.3) 
and all other skeleton loop integrals are given explicitly below as follows: 

'j^T'iiq, ^qi,qi - q)T3{-q, qi,q - qi)D{qi)D{q ~ gi), (2.4) 
T2{q) ^ g'^ J J ^^^;^T^T3{~q2,q3,q2 - q3)T3i~q,qi,q3 ~ q2)D{qi)D{-q2)D{q3)D{q3 ~ q2), (2.5) 

niq)=9'l ^ / "^nU~q,q^,^q2,q3)D{q^)Di-q2)Diq3), (2.6) 

where in the last two skeleton loop integrals q — qi + q2 — q3 — as usual. Evidently, neither the color group factors 
nor the Dirac indices play any role in tracking down the mass gap, which can only be of dynamical origin. That is 
the reason why they are not explicitly shown throughout this paper. 

The general iteration solution (i.e., when the above-mentioned skeleton loop integrals are to be iterated) of the 
gluon SD equation (2.1) looks like 

D{q) = D^°'>{q) + D^^'>{q) + D^^'>{q)+D^^\q) + ... 

= D"{q) + D°{q)Tg[D'> + D^'^ + + ^(3) + ,„](q)[D''{q) + D^^'> (q) + D^^\q) + D^^\q) + ...], 

(2.7) 

and £)'"'((/) = D^{q)- It is nothing but the skeleton loops expansion, and it is not the PT series. First of all, the 
magnitude of the coupling constant squared cannot be fixed to be small, i.e., it is arbitrary. Secondly, the dependence 
of the skeleton loop integrals on the coupling constant squared is also completely arbitrary, i.e., it cannot be explicitly 
fixed on general ground. However, this expansion is rather formal, since the corresponding skeleton loop integrals are 
not yet regularized. In the deep IR limit q^ ^ the skeleton loop integrals (2.4)-(2.6) tend to their corresponding 
divergent constant values (in Euclidean metrics g^ — > implies g^ — > 0). Just these constants having the dimensions 
of a mass squared are the main objects of our investigation. The only problem is as how to extract them, i.e., to 
make their existence and important role perfectly clear. Fortunately, this can be explicitly done through the necessary 
regularization of the initial skeleton loop integrals. 
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III. REGULARIZATION 

Due to AF P, 01 all the skeleton loop integrals as well as those which will appear in the formal iteration solution 
(2.7) are divergent, so the general problem of their regularization arises. For our future purpose it is convenient to 
regularize the above-mentioned skeleton loop integrals by subtracting as usual their values at a safe (slightly different 
from zero) space-like point = jj? (Euclidean signature is already chosen). Thus, one obtains 

T^{q) = T,{q)-T,{q^^^,''), 
T^{q) = T2{q)-T2{q^^^i''), 

= n{q)-n{q^^^^^), (3.1) 

where all the divergent constants having the dimensions of a mass squared are 

Ti{q^ = H^) = j ^^^^T^{q,-qi,qi - q)T3{-q,qi,q - qi)D{qi)D{q - qi)\g2^^2, 

I 1 i\ 4 f id*qi f i(Pq2 o / 
T2(q ^^i)^g J j—j^ J j^;^T^T3{-q2,q3,q2~q3)T3{-q,qi,q3~q2) 



D{qi)D{-q2)D{q3)D{q3 - q2)\,2^_ 



T^{q' ^ ,,^) ^ I ^ J '-^^T2U-q,q,,-q2,q3)D{q,)D{~q2)D{q3)\,2^,2. (3.2) 

At the same time, the introduction of the ultraviolet (UV) cutoff A^ in these integrals is assumed . All this makes it 
possible first of all to explicitly release the above-mentioned divergent constants, and put all the divergent skeleton 
loop integrals under firm control. The UV cutoff should go to infinity at the final stage only. The formal subtraction 
procedure for the constant tadpole term (2.3) implies = Tt — Tt ^ 0. This is in agreement with the dimensional 
regularization method HQ where Tt{D°) = 0, indeed. So, we can discard this term within the formal iteration 
solution (2.7) without loosing generality. 

The subtractions in Eq. (3.1) mean that the decomposition of the regularized quantities into the independent tensor 
structures can be written down as follows: 

T«(g) ^ T« ^,(q) ^ 6,.q'T^'\q') + q,qM'\<l% 
T?iq) ^ T« ^,(q) = 5,,q^Ti'\q^) + q,qM^\q% 
n"i<i) ^ Tl^),M = 6,.q'T^^'\q') + q,q.T^^'\q'). (3.3) 

where all invariant functions T^'"\q'^),T^''\q'^),T^''"\q'^) at n ^ 1,2 are dimensionless ones. In the region of small 
q^ they are represented in the form of the corresponding Taylor expansions and remain arbitrary otherwise. Due to 
the definition 5^5,^ — q^L^^, instead of the independent structures S^i, and q^qv in Eqs. (3.3) and below, one can 
use and L^^^ as the independent structures with their own invariant functions. From these relations it follows 
that the regularized skeleton loop contributions are always of the order g^, i.e., T^{q) = 0{q^), T^{q) — 0{q^) and 
n^iq) = 0(q2). 

The NL pure gluon part (2.2) thus can be exactly decomposed as the sum of the two terms: 



where 



and 



T,mq)^T,[D]+T^^[D]{q), (3.4) 



T,[D] = ^UD') + lT2iD^) + IniD'), (3.5) 



T.'^'mq) - iTf iq) + ^T2%) + (g) = 0{q';D). (3.6) 
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In Eq. (3.5) we introduce the following notations: Ti{q^ = y?) = Ti{D^), T2[q^ = n^) = T2{D^) and T^(g2 ^ = 
T2{D^), showing their dependence on the corresponding number of the gluon propagators only. In Eq. (3.6) the last 
equality shows explicitly that its left-hand-side is always of the order q'^, depending again on D in general. Thus the 
gluon SD equation (2.1) becomes divided into the two principally different terms, namely 

D{q) = D'^iq) + D\q)T,[D]D{q) + D\q)0{q^ ■ D)D{q). (3.7) 

Let us note that both blocks, Tg[Z)] and 0{q^;D) (which should be iterated with respect to D), depend on D 
themselves. This means that there is no way to sum up in general these iteration series into the corresponding 
geometric progression. The latter possibility appears only when the blocks to be iterated do not depend on D like in 
quantum electrodynamics (QED), where the electron skeleton loop is only present. The difference between these two 
cases lies, of course, in the NL dynamics of QCD in comparison with the linear one in QED. 

Let us emphasize once more that the constant block Tg[£'] having the dimensions of a mass squared is just the 
object we have worried about to demonstrate explicitly its crucial role within our approach. In this connection a few 
additional remarks are in order. As it follows from the standard gluon SD equation (3.7), the corresponding equation 
for the gluon self-energy looks like 

D-\q) = D^\q) - Tg[D] - 0{q'-D), (3.8) 

where we put D^{q) = Do(g). In order to unravel overlapping UV divergence problems in YM theory, the neces- 
sary number of the differentiation with respect to the external momentum should be done first (in order to lower 
divergences). Then the point-like vertices, which are present in the skeleton loop integrals, should be replaced by 
their full counterparts via the corresponding integral equations. Finally, one obtains the corresponding SD equations 
which are much more complicated that the standard ones, containing different scattering amplitudes, which skeleton 
expansions are, however, free from the above-mentioned overlapping divergences. Of course, the real procedure [Tll | 
(and references therein) is much more tedious than briefly described above. However, even at this level, it is clear 
that by taking derivatives with respect to the external momentum q in Eq. (3.8), the main initial information due 
to the constant block will be totally lost. Whether it will be restored somehow or not at the later stages of 

the renormalization program is not clear at all. Thus in order to remove overlapping UV divergences ("the water") 
from the SD equations and skeleton expansions, we are in danger to completely loose the information on the constant 
block Tg[D] which is the dynamical source of the mass gap ("the baby") within our approach. In order to avoid 
this danger and to be guaranteed that no any dynamical information are lost, we are just using the standard gluon 
SD equation (3.7). The presence of any kind of UV divergences (overlapping and usual (overall)) in the skeleton 
expansions will not cause any problems in order to detect the mass gap responsible for the IR structure of the true 
QCD vacuum. In other words, the direct iteration solution of the standard gluon SD equation (3.7), complemented 
by the proposed regularization scheme, is reliable to release a mass gap, and thus to make its existence perfectly clear. 
The problem of convergence of such regularized skeleton series is completely irrelevant in the context of the present 
investigation. Anyway, we keep any kind of UV divergences under control within our method. At the same time, any 
kind of UV divergences play no any role in the existence of a mass gap responsible for the IR structure of the full 
gluon propagator, i.e., its existence does not depend on whether overlapping divergences are present or not in the SD 
equations and corresponding skeleton expansions. All this is the main reason why our starting point is the standard 
gluon SD equation (2.1) for the unrcnormalized Green's functions (this also simplifies notations). 

IV. NONLINEAR ITERATION 

The formal iteration solution of the gluon SD equation (3.7) now looks like 

D{q) D^''\q) + D^^^{q) + D^^Hq) + D^''^{q) + ... = D''{q) 

+ D''iq)Tg[D'' + D^^'> + + ^(3)(g) + ...][D°{q) + D^^'>{q) + D^^^q) + D^^'Hq) + ...] 
+ D''iq)0{q^;D° + D^^'> + D^^'> + D^^\q) + ...)[D°{q) + D^^\q) + D^''\q) + D^^\q) + ...], 

(4.1) 



where 
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- D"iq)F,[D'>]D"{q) + D%q)0{q';D'>)D'{q), 

D^^\q) = D"{q)F2[D'>,D^^'>][D°iq) + D^^'>{q)]+D"{q)0{q^;D° + D^'^)[D°{q)+D^'\q)], 

D(^\q) = D"{q)F^[D°,D^^\D^^^]lD'>{q)+D^^\q) + D^^'>{q)] 

+ D"{q)0{q^;D° + D^'^ +D'^^y)[D"{q) + D^^'>{q)+D^^\q)], (4.2) 

and so on. Let us consider the divergent constants Fi[D^], F2[D^ , D^^^], F^ID^ , D^^^ JD^^^], introducing the following 
notations: 

F,[D°] = T,[D°] = ^T,i{Dy) + ^n{{Dy) + ^n{{Dy), 

F2[i?",i?(i)] = TgiDO + DW] = lri((i?" + i?(i))2) + iT2((i?o + L»«)4) + ir^((Z)0 + Z)(i))3), 

2 6 

and so on. As underlined above, each of them has the dimensions of a mass squared, so on general ground one can 
represent them as follows: 

F,=F4D°] - A'C,{X,u,^,g'), 

F2^F2[D\D^^^] = A2C2(A,j.,C,5'), 

F3^i^3[i?°,^('\i?('^] = A2C3(A,z.,C,5'), (4.4) 

and so on. In these relations is the above-mentioned mass gap. The dimensionless constants Ci,C2,C3,... 
depend on the dimensionless UV cutoff A introduced as follows: = AA^. They depend on the dimensionless 

renormalization point v introduced as follows: /i^ = z^A^. Evidently, via the corresponding subscripts these constants 
depend on which iteration for the gluon propagator D is actually done, which in its turn gives the dependence on the 
gauge-fixing parameter ^. They also depend on the dimensionless coupling constant squared (see Eqs. (3.2)). The 
parameters A and ;/ may depend on ^ and (and hence vise- versa). The dependence of the mass gap A^ on all these 
parameters is not shown explicitly, for convenience, but can be restored any time, if necessary. 

In the relations (4.3) and (4.4) we use the short-hand notation = D^{q), however, it is more appropriate to 
introduce the short-hand notations as follows: 



D\q) ^ Doiq)^Do, 0,{q^) ^ 0{q^-,D\ 

02{q^) = 0{q^; D° + D^'^), 0^{q^) = 0{q''; + D^'^ + D^^)), (4.5) 
and so on. Then the formal iteration solution (4.1) becomes 

D{q) = DW(g) + D^^\q) + D^^^q) + D^^\q) + ... = Do + [DoF^Dq + DoOi{q^)Do] 

+ [D0F2D0 + DaF2DoFiDo + DoF2DoOi{q'')Do 

+ DMq')Da + Do02{q^)DoFiDo + DMq^)DoOi{q^)Do] 

+ [DoFgi^o + DoF:,DoFiDa + DaF^D,,Oi{q^)Do 

+ D0F3D0F2D0 + D0F3D0F2D0F1D0 + DoF3DoF2DoOiiq^)Do 

+ DoF3Do02{q'^)Do + DoF3Do02{q'^)DoFiDo + DoF3Do02iq^)DoOi{q'^)Do 

+ Do03iq'')Do + Do03iq'')DoF,Do + DMq^)DM<l^)Do 

+ Da03{q^)DaF2Do + Da03iq^)DaF2DoFiDo + Do03{q^)DoF2DoOi{q^)Do 

+ Do03{q^)Do02{q^)Do + Do03{q^)Do02{q^)DoF^Do + Da03{q^)Do02{q^)DQOi{q^)Do] + ... . 

(4.6) 
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V. SHIFTING PROCEDURE 

The formal iteratfon series (4.6), however, much more convenient to equivalently rewrite as follows: 

D{q) = [Dl{Fi +F2+F3 + ...) + Dl{FiF2 + FsFi + F2F3 + ...) + D^iF^FiF^ + ...) + ...] 

+ [DoF2DoOi{q'')Do + Do02{q'')DoFiDo + DoF3DoOi{q^)Do + DoF3DoF2DoOi{q^)Do 

+ DoF3Do02{q^)Da + DoF^Do02{q^)DoFiDo + DoF3Do02{q^)DoOi{q^)Do 

+ Do03iq^)DoFiDo + Do03{q^)DoF2Do + Do03{q^)DoF2DoFiDo 

+ Do03{q'')DoF2DoOiiq'')Do + Do03{q^)Do02iq^)DoFiDo + ...] 

+ [Do + DoOi{q^)Do + Do02{q'')Do + Do03{q')Do + Do02{q^)DoO^{q'')Do 

+ Da03{q^)DoO^{q^)Da + Do03{q^)Do02{q^)Do + Do03{q^)Da02{q^)DoOr{q^)Do + ...]. 

(5.1) 

Since Dq = Do{q) there arc three formally different types of terms. The first type of terms is singular 

as much as D^'^'^ ~ (Q'^)~^~'^, k = 0,1,2,3,.... Evidently, only the divergent constants Fi,F2,F3,... and their 
combinations enter into these terms. The third type of terms is singular as much as the free gluon propagator, i.e., 
they arc of the order Do ^ (9^)^^) since all functions On{q'^),n = 1,2,3, ... are of the order q'^. Evidently, only the 
functions On{q^),n — 1,2,3, ... and their products enter into these terms. The second type of terms is the so-called 
mixed up terms which contain the divergent constants ^1,^2,^3,... and the functions 0„(g^),n = 1,2,3,... in the 
different combinations. 

Let us show that any mixed up term is simply the exact sum of the first and third types of terms. In order to show 
this explicitly, let us recall that all functions 0„((jr^), n = 1, 2, 3, ... arc regular functions of q^ . So, wc can equivalently 
represent them in the form of the corresponding Taylor expansions. It is convenient to present such kind of expansions 
in terms of D^^ ~ q^ and making an exact decomposition as follows: 

fe-1 fe-1 

On(g') = O„(a')-^0"'E(^0^')""^^^(l) + ^0'E(^0^')"'"On'"Hl) 

m=0 m=0 

k-1 

= Do-iE(^o^')""Oi"Hl) + Doi(DoA2)-V„(g2)^ n= 1,2,3,... . (5.2) 

TO=0 

Here fn{Q^) are dimensionless functions with constant behavior at zero momentum and otherwise remaining arbitrary. 
The last term always has the corresponding order in powers of Dq^, so that 

fe-1 

D^'iDoAY'fniq') = 0„(g^)-Do"'E(^o^')""^"'"^W 

m=0 

00 fc— 1 

= Do"' E (^oA^)— Oi'")(l) - Dq-i E (^oA2)-"Oi™)(l) ~ Do\DoAY', 

m=0 m=0 

(5.3) 

indeed. The same is true for the products of two, three, and more different On{q^) functions (see below). 
For the first mixed up term, DoF2DQOi{q^)Do, one then gets 

DoF2DoO,{q^)Do = 0^+' F2[D^'o'f\l) + D^\DoA')-' f,{q^)] 

= DlF20f\l) + DoF2A-^h{q^) = DlF20f\l) + 0{Do), (5.4) 

since for this term we should put fc = 1 in Eq. (5.2). The combination of the mass squared parameters _F2A^^ is the 
dimensionless one (see relations (4.4)), and the last term is obviously denoted as 0{Dq). Thus, one concludes that 
this mixed up term becomes the exact sum of the two different terms. The first one should be included into the terms 
of the order Dq shown as the first term in the iteration solution (5.1), while the second term should be combined with 
the terms which always are of the order Dq. 
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For the mixed up term, Do^3-Do^2-DoOi (9^)-Do, we obtain 

DoF,DoF2DoOiiq^)Do = DI+^ F^F^lD^' of\l) + D^^^-^O^l\l) + D^HD,^^)-^ h{q^)] 

= DlF2F:,0f\l) + DlF2F:,^-^0^^\l) + DaF^F2^-^ hiq") 

= DlF2F:,0f\l) + DlF2F:,l:,-^0't\l) + Opo), (5.5) 

since for this term we should put A: = 2 in Eq. (5.2). The combination F2F^/S.^^ has the dimensions of a mass squared 
(see Eqs. (4.4)), while the combination F2F^/S.''^ is dimensionless. Again the last term is denoted as 0(D"). So this 
term also becomes the exact sum of the three terms. The first term from this sum has to be shifted into the second 
term of Eq. (5.1). The second term from this sum has to be shifted into the first term of Eq. (5.1), and the last one 
has to be shifted into the last term of the general iteration solution (5.1). 

It is instructive to consider in more details the terms which contain two and more 0{q^) functions. The corresponding 
Taylor expansions for the product of any two and three functions are as follows: 

m=0 
fc-3 

Or.{q^)OM^)0,{q^) = D^'Y.iDoAY"'oi:f{l) + D,'{DoA')-'^'fniA<l'), (5-6) 

m=0 

and so on. Here fni{q^) and fnij{q^) are dimensionless functions with constant behavior at small momentum q^ and 
otherwise remaining arbitrary. Similarly to the previous case, the last terms are the terms of the corresponding orders 
in powers of D^^^ so that 

I?o-2(i^oA2)-'=+i/„K9') - 0„K^o'"')-^o'PoA2)-'=+i, 

D^\DoA')-^+^ UiA^^) ^ Oru,{D^'-')-D^'{DoA')-''+', (5.7) 
indeed, and so on. Then, for example, for the mixed up term, Do03{q^)DoF2DoOi{q^)Do, one gets 

Do03{q^)DoF2DoOi{q^)Do = Dl+^F2[D^^O^°,\l) + D^\DoA^)-^f3,{q^)] 

= DlF20^°,\l)+DoF2A-^fsi(/) 

= DlF20^^^{l) + OiDo), (5.8) 

since for this term we should put /c = 2 in the first of Eqs. (5.6). The combination F2A~^ is dimensionless, and the 
last term is denoted as 0{Dq). Again this mixed up term becomes the exact sum of two terms. The first term should 
be shifted into the first, while the second one should be shifted into the last term of the general iteration solution 
(5.1). 

Moreover, all other mixed up terms, which are explicitly present and omitted in formal series (5.1), should be 
treated in the same way. Let us underline that the exact and unique separation between the two kind of terms 
(~ D^'^^{q)^ fc = 0, 1, 2, 3, ... and ~ Do{q)) is achieved by keeping the necessary number of terms in the corresponding 
Taylor expansions. This "shifting" method in its general form has been formulated and applied in our previous 
publications 0, ^3 . It is worth emphasizing that by shifting we do not change the functional dependence (and hence 
the dependence on the mass gap) in the terms ~ Z)g^'^(q), only the accompanied q^-independent factors will be 
changed. At the same time, the terms ~ DQ{q)) will be changed (many new terms of the same order will appear). 
Let us also note that all the combinations of different masses squared (like the above-mentioned F2A^^, F2F3A^^, 
etc.) can be reduced to the corresponding powers of the mass gap A^ and dimensionless coefficients Ci, C2, C3, ... via 
the relations (4.4). They are to be multiplied by the different dimensionless constants On"'\l), ol^\l), Ol^j{l), 
etc., which appear through the shifting procedure, and which in general may depend on the same quantities as in Eq. 
(4.4). 

Rearranging all the terms, one gets that the general iteration solution (5.1) for full gluon propagator becomes the 
exact sum of the two different terms, namely 

00 00 

D{q) - Dl{q)A^ Y.^D,{q)A^]>^ ^ a,..„(A, i., e, g') + 0{Do{q)). (5.9) 

fc=0 m=0 
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It is worth emphasizing once more that Eq. (5.9) obtained by the shifting procedure is equivalent to the initial 
general iteration series (5.1). In other words, the shifting method is not an approximation, but the exact method 
of the corresponding rearrangement of the terms. It makes it possible to represent initial series (5.1) in the form 
much more convenient for our purpose (to track down the mass gap), while shifting functional ambiguity of the initial 
series (5.1) (which is due to the functions 0„(g^)) into the term 0{Do{q)) in Eq. (5.9). Let us recall that the terms 
0{Do{q)) is the sum of the terms which are of the order DQ{q) from the very beginning (the third type of terms in 
the expansion (5.1)) and the terms of the same order having appeared due to the above-described shifting procedure. 

VI. EXACT STRUCTURE OF THE FULL GLUON PROPAGATOR 

Restoring the tensor structure, omitting the tedious algebra and again taking into account that Do{q) ^ {q'^)~^, 
the general iteration solution of the gluon SD equation (5.9) for the full gluon propagator can be algebraically (i.e., 
exactly) decomposed as the sum of the two principally different terms as follows: 

A 2 oo oo oo -J 



(q2)2 



where the superscript "INP" stands for the intrinsically NP part of the full gluon propagator. Its exact structure 
inevitably stems from the general iteration solution of the standard gluon SD equation. The important feature of our 
method is that the skeleton loop integrals have been iterated (skeleton loops expansion), so no any assumptions and 
approximations have been made. We distinguish between two terms in Eq. (6.1) by the character of the corresponding 
IR singularities and the explicit presence of the mass gap (see below). It is worth emphasizing that both terms are 
valid in the whole energy/momentum range, i.e., they are not asymptotics. At the same time, we achieved the exact 
separation between two terms responsible for the NP (dominating in the IR (g^ 0)) and PT (dominating in the 
UV {q^ — > oo)) dynamics in the true YM vacuum. 

The PT part of the full gluon propagator remains undetermined. The exact dependence of the PT gluon form factor 
dF^{q^,£) = J2m=o '^"lil'^ ' oil 9^ cannot be fixed on general ground like it has been done in its INP counterpart 
(what we only know about 0^(9^, C)"functions is that all of them are regular at small q^, and the sum over them 
produces the PT logarithm improvements at large q^ due to AF). Evidently, the presence of overlapping and overall 
UV divergences in the PT part of the full gluon propagator cannot change the structure of its INP part, i.e., its 
functional dependence on and hence its dependence on the mass gap. They may only affect the g^-independent 
factors 0fc^m(A, J^, ^, g^), which concrete values, however, are not important. In the PT part the sum over m indicates 
that all iterations contribute into the PT IR singularity only, which is defined as always being as much singular as 
the power-type IR singularity of the free gluon propagator {q^)~^. That is why the longitudinal component of the full 
gluon propagator should be included into its PT part. Anyway, we are not responsible for this part. It is the 
have payed to fix exactly the functional dependence of the INP part of the full gluon propagator. In Refs. 
we came to the same structure (6.1) but in a rather different ways. 

A. The INP phase in QCD 

The exact decomposition of the full gluon propagator into the two principally different terms in Eq. (6.1) is only 
possible on the basis of the corresponding decomposition of the full gluon form factor d{q'^,£^) in Eq. (1.1), namely 

d(9^c) = d(<z^o-rf^^('?^e)+d^^(<z^^) = rf"^^(9^c)+rf^^(<z^o■ (6.2) 

Let us note that in principle the full gluon form factor can be defined as the effective charge of QCD, i.e., d{q^) = as{q'^), 
where the dependence on ^ is omitted, for simplicity. This algebraic decomposition makes it possible to define correctly 
the NP phase in comparison with the PT one. The full gluon form factor d(q'^,£,) being the NP effective charge, 
nevertheless, is "contaminated" by the PT contributions, while d^^{q^, ^) is the truly NP one, since it is free of them, 
by construction |j7|,|lfl|. Substituting the exact decomposition (6.2) into the full gluon propagator (1.1), one obtains 



DM ^^{TMd{q^O+^LM} 1 = Dj^^^iq)+Dj:^{q), (6.3) 
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where 



and 



q" 



(6.5) 



in complete agreement with Eq. (6.1), indeed. 

As it follows from Eq. (6.1), the INP part of the full gluon propagator in Eq. (6.4) is nothing else, but the 
corresponding Laurent expansion in integer powers of accompanied by the corresponding powers of the mass gap 
squared and multiplied by the sum over the ^^-independent factors, namely 



oc 

i^^r (9, A2) = ^T,M J2{qT'-\^')'^'MK ^, 5^), (6.6) 

k=0 

where 



M\ 9^) = E '^fc^"(^' ^' 5') (6.7) 

are the so-called residues at poles. The sum over m indicates that an infinite number of iterations (all iterations) of 
the corresponding skeleton loop integrals invokes each severe (for definition see below) IR singularity labelled by k. 

Using the exact decomposition of the full gluon propagator described above, we can define in general terms the INP 
phase in QCD as follows: 

(i) . The INP phase is characterized by the presence of the power-type severe (or equivalently NP) IR singularities 
(9^)^^ '^i k — 2, 3, .... So these IR singularities are defined as more singular than the power-type IR singularity of 
the free gluon propagator (g^)^^. The Laurent expansion (6.6) necessarily starts from the simplest NP IR singularity 
(g^)~^ possible in four-dimensional QCD, indeed y\. 

(ii) . It depends only on the transverse ("physical") degrees of freedom of gauge bosons. 

(iii) . It is gauge- invariant. Though the coefficients (f)k,rni^,i^,£,i 9^) of the Laurent expansion (6.6) may explicitly 
depend on the gauge-fixing parameter ^, the structure of this expansion itself does not depend on it (see discussion 
below as well). 

(iv) . The INP part of the full gluon propagator vanishes as the mass gap goes to zero, while the PT part survives. 

Within our approach the mass gap determines the power-type deviation of the full gluon propagator from the free 
one in the IR limit (g^ 0) , while Aqc d determines the logarithmic deviation of the full gluon propagator from the 
free one in the UV limit (g^ — > 00). So, we distinguish between the two different phases in QCD not by the strength 
of the coupling constant squared (which is arbitrary in our approach), but rather by the explicit presence of the mass 
gap, in which case the coupling constant plays no any role. The INP phase disappears when it goes to zero even if the 
NP IR singularities are not explicitly present. So the subtraction (6.2) can be equivalently written down as follows: 



d^^(g2, A^) = d(g2, A^) - d(g2, A^ = 0) = d(g2, A^) - d^^iq^), (6.8) 

i.e., it goes to zero as A^ —> 0, indeed. This once more emphasizes the important role of the mass gap in the definition 
of the truly NP phase as a particular case of the INP one when the NP IR singularities are not explicitly present as 
stated above. In other words, the existence of the truly NP phase in any approach based on the gluon propagator 
assumes the regular dependence on the mass scale parameter, which is chosen to play the role of the mass gap. 
Otherwise, in the absence of the mass gap and in order to recover the truly NP phase the UV asymptotic of the full 
gluon propagator should be subtracted in agreement with the relation (6.8). It is worth emphasizing once more that 
the existence of the mass gap and the presence of the NP IR singularities in the Laurent expansion (6.6) (and hence 
in the full gluon propagator (6.1)) is absolutely general phenomenon. It does not depend on the concrete values of 
the parameters: A, i/, ^, g^. It is only due to the NL interaction of massless gluons. 
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VII. DISCUSSION 

The unavoidable presence of the first term in Eq. (6.1) makes the principal difference between non-abelian QCD 
and abelian QED, where such kind of term in the full photon propagator is certainly absent (in the former theory 
there is direct coupling between massless gluons, while in the latter one there is no direct coupling between massless 
photons). Precisely this term may violet the cluster properties of the Wightman functions [ij, and thus validates 
the Strocchi theorem which allows for such IR singular behavior of the full gluon propagator. 

The INP part of the full gluon propagator in the form of the corresponding Laurent expansion describes the so-called 
zero momentum modes enhancement (ZMME or simply ZME which means zero momentum enhancement) effect in 
the true QCD vacuum due to the NL dynamics of massless gluon modes there. As underlined above, we do not specify 
explicitly the value of the gauge-fixing parameter f . So the ZMME effect takes place at its any value. In this sense 
this effect is gauge-invariant. This is very similar to AF. It is well known that the exponent which determines the 
logarithmic deviation of the full gluon propagator from the free one in the UV region (g^ ^ ^%cd) explicitly depends 
on the gauge-fixing parameter. At the same time, AF itself does not depend on it, i.e., it takes place at any ^. 

Also, due to the arbitrariness of the above-mentioned residues 0fc(A, v, ^, g^) there is no way to sum up these Laurent 
series into the function with regular behavior at small . However, the smooth in the IR gluon propagator is also 
possible depending on different truncations/approximations used, since the gluon SD equation is highly nonlinear one. 
The number of solutions for such kind of systems is not fixed a priori. The singular and smooth in the IR solutions 
for the gluon propagator are independent from each other, and thus should be considered on equal footing. Anyway, 
in order to find the smooth gluon propa gat or completely different (from direct iteration) method of the solution of 
the gluon SD equation should be used [8|, Il4| (and references therein as well) , since we have explicitly shown here in a 
gauge-invariant way and making no any approximations/truncations that the general iteration solution is inevitably 
severely singular at small gluon momentum. 

The QCD Lagrangian does not contain a mass gap. However, we discovered that the mass scale parameter re- 
sponsible for the NP dynamics in the IR region exists in the true QCD ground state. At the level of the gluon SD 
equation it is hidden in the skeleton loop contributions into the gluon self-energy. It explicitly shows up (and hence 
the corresponding severe IR singularities) when the gluon momentum goes to zero. At the fundamental quark-gluon 
(i.e., Lagrangian) level the dynamical source of the mass gap is the triple and quartic gluon vertices, i.e., the NL 
dynamics of QCD. The former vanishes when all the gluon momenta involved go to zero (r3(0,0) = 0), while the 
latter survives in the same limit (T4(0,0, 0) ^ 0). Because of these features it would be tempting to think that the 
quartic potential {A A A)^ in the action plays much more important role than its triple counterpart in the arising of 
severe IR singularities in quantum YM theory. However, since we are dealing with the skeleton loops expansion, we 
are unable (at least at this stage )to exactly establish that the mass gap could arise from the quartic potential rather 
than from its triple counterpart [a] (see also Ref. 15]). So this dilemma remains an open but an interesting problem 
to be solved. 

In this connection, it is necessary to discuss an interesting feature of the INP part of the full gluon propagator. Its 
functional dependence on and hence on the mass gap (i.e., the Laurent structure of the expansion in Eq. (6.6)) 
is exactly established up to the corresponding residues (/>fe(A, z/, ^, g^). In these residues the contributions from both 
the triple and quartic gluon vertices have been taken into account. However, the Laurent structure of the INP part 
does not depend on whether we will take into account only three- or only four-gluon vertices or both of them in 
the above-mentioned skeleton loop integrals. We cannot omit both, so one of the NL interactions should be always 
present. On the other hand, the residues themselves will depend, of course, on the character of the NL interaction 
taken into account. However, obviously, their concrete values are not important, only the general dependence of the 
residues (?I>fc(A, J^, Ci5^) o^i their arguments is all that matters (in the subsequent paper we will show this explicitly). 
At the same time, the dependence of the functions 0^(9^; ^) and hence of the PT form factor d^'^{q^; ^) in Eq. (6.1) 
on will heavily depend on the character of the NL interactions taken into account in the corresponding skeleton 
loop integrals but this is not important for us as underlined above. 

Thus the true QCD vacuum is really beset with severe IR singularities. They should be summarized (accumulated) 
into the full gluon propagator and effectively correctly described by its structure in the deep IR domain, exactly 
represented by its INP part. The second step is to assign a mathematical meaning to the integrals, where such kind 
of severe IR singularities will explicitly appear, i.e., to define them correctly in the IR region 0,01 ^3- This can be 
done by the use of the dimensional regularization method correctly implemented into the distribution theory [T^ 
(see subsequent paper). Just this IR violent behavior makes QCD as a whole an IR unstable theory, and therefore it 
has no IR stable fixed point, indeed | l|. [l p . This means that QCD itself might be a confining theory without involving 
some extra degrees of freedom [H llSl Il9| fand references therein). 

There is no doubt that our solution for the full gluon propagator, obtained at the expense of remaining unknown its 
PT part, nevertheless, satisfies the gluon SD equation (3.7), since it has been obtained by the direct iteration solution 
of this equation. To show this explicitly by substituting it back into the initial gluon SD equation is not a simple 
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task, however, and this is to be done elsewhere. The problem is that the decomposition of the full gluon propagator 
into the INP and PT parts by regrouping the so-called mixed up terms in section V was a well defined procedure 
(there was an exact criterion introduced in section VI how to distinguish between these two terms in a single D). 
However, to do the same at the level of the gluon SD equation itself, which is nonlinear in D, is not so obvious. Also, 
the corresponding severe IR singularities should be put under control at first within the distribution theory (all this 
will be explicitly demonstrated in a forthcoming paper). 

VIII. CONCLUSIONS 

A few years ago Jaffe and Witten (JW) have formulated the following theorem P|: 

Yang-Mills Existence And Mass Gap: Prove that for any compact simple gauge group G, quantum Yang-Mills 
theory on exists and has a mass gap A > 0. 

Of course, at present to prove the existence of the YM theory with compact simple gauge group G is a formidable 
task yet. It is rather mathematical than physical problem. However, one of the main results of our investigation here 
can be formulated similar to the above-mentioned JW theorem as follows: 

Mass Gap Existence: If quantum Yang-Mills theory with compact simple gauge group G = SU{3) exists on R^, 
then it exhibits a mass gap. 

Our mass gap remains neither IR nor UV renormalized yet, since at this stage it has been only regularized, i.e., 
A^ = A^(A, ^, g^). However, there is no doubt that it will survive both renormalization programs (see subsequent 
paper). So denoting its IR and UV renormalized version in advance as A^vp, then a symbolic relation between it, the 
JW mass gap (A = Ajw) and Aqcd = Apt could be written as 

Anp < — (u^MiR M^v^oo — * ^PT- (8.1) 

Here as is obviously the fine structure coupling constant of strong interactions, while Mjjv and Mjn are the UV and 
IR cut-offs, respectively. The right-hand-side limit is well known as the weak coupling regime, while the left-hand-side 
can be regarded as the strong coupling regime. We know how to take the former 1, 2, 18], and we hope that we have 
explained here how to deal with the latter one, not solving the gluon SD equation directly, which is a formidable task, 
anyway. However, there is no doubt that the final goal of this limit, namely, the mass gap A^p exists, and should be 
renormalization group invariant in the same way as Aqcd- It is solely responsible for the large scale structure of the 
true QCD ground state, while Apt is responsible for the nontrivial PT dynamics there. 

It is important to emphasize once more that the mass gap has not been introduced by hand. We have explicitly 
demonstrated that it is hidden in the skeleton loop integrals contributing into the gluon self-energy due to the NL 
interaction of massless gluon modes (Eqs. (2.4)-(2.6)). The mass gap shows explicitly up when the gluon momentum 
goes to zero. An appropriate regularization scheme has been applied to make the existence of the mass gap perfectly 
clear. Moreover, it survives an infinite series summation of the corresponding skeleton loop contributions (skeleton 
loop expansion). In other words, an infinite number of iterations of the relevant skeleton loops has to be made in order 
to exhibit a mass gap. No any truncations/approximations have been made as well as no special gauge choice, i.e., the 
result of the summation is gauge-invariant. In the presence of the mass gap the QCD coupling constant plays no any 
role. All its orders contribute into the mass gap (skeleton loop expansion). This explains why the interaction in our 
picture can be considered as a strong one. It is worth emphasizing that our mass gap and the JW mass gap cannot 
be interpreted as the gluon mass, i.e., they always remain massless. These features point out on a possible similarity 
between our mass gap and the JW one. Moreover, in a next paper we will explicitly show (by investigating the IR 
renormalization properties of the mass gap) that the interaction in our picture is not only strong but short-ranged as 
well. This will allow us to analytically formulate the gluon confinement criterion in a gauge-invariant way. 

Our second important result, a byproduct of the proof of the existence of a mass gap, is that the general iteration 
solution of the gluon SD equation is unavoidably severely IR singular at small gluon momentum. The exactly 
established Laurent structure of the INP part of the full gluon propagator (6.1) clearly shows that an infinite number 
of iterations of the relevant skeleton loops (skeleton loops expansion) invokes each NP IR singularity. Again no special 
gauge choice and no any truncations/approximations have been made as well in such obtained general iteration solution 
of the gluon SD equation for the full gluon propagator. So, our gluon propagator (more precisely its INP part) takes 
into account the importance of the quantum excitations of severely singular IR degrees of freedom in the true QCD 
vacuum. They lead to the formation of the purely transverse quantum virtual field configurations with the enhanced 
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low-frequency components/large scale amplitudes due to the NL dynamics of the massless gluon modes. We will call 
them the purely transverse singular gluon fields, for simplicity. In the presence of such severe IR singularities the IR 
multiplicative renormalization (IRMR) program is needed to perform in order to remove them in a self-consistent way 
from the theory (see subsequent paper). In another forthcoming paper we will show that the quark and ghost degrees 
of freedom play no any significant role in the dynamical generation of a mass gap. The NL interaction of massless 
gluon modes is only important within our approach. 

In summary, the existence of a mass gap in quantum YM theory has been proved in a gauge-invariant way. We 
have explicitly shown (again in a gauge-invariant way) that the direct iteration solution of the gluon SD equation 
is inevitably severely singular at small gluon momentum. Our general conclusion is that the behavior of QCD at 
large distances is governed by the mass gap, and therefore it should play a crucial role in the NL realization of the 
quantum-dynamical mechanism of color confinement. 
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